Semilocal and Hybrid Density Embedding Calculations of Ground-State 

Charge-Transfer Complexes 
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We apply the frozen density embedding method, using a full relaxation of embedded densities 
through a freeze-and-thaw procedure, to study the electronic structure of several benchmark ground- 
state charge-transfer complexes, in order to assess the merits and limitations of the approach for 
this class of systems. The calculations are performed using both semilocal and hybrid exchange- 
correlation (XC) functionals. The results show that embedding calculations using semilocal XC 
functionals yield rather large deviations with respect to the corresponding supermolecular calcula- 
tions. Due to a large error cancellation effect, however, they can often provide a relatively good 
description of the electronic structure of charge-transfer complexes, in contrast to supermolecular 
calculations performed at the same level of theory. On the contrary, when hybrid XC functionals 
are employed, both embedding and supermolecular calculations agree very well with each other and 
with the reference benchmark results. 

In conclusion, for the study of ground-state charge-transfer complexes via embedding calcula- 
tions hybrid XC functionals are the method of choice due to their higher reliability and superior 
performance. 

PACS numbers: 



INTRODUCTION 



(KSCED) [2J,l2£ 



In the last years a considerable interest in density- 
embedding methods has led to a widespread development 
of subsystem approaches |ll421 1 within density functional 
theory (DFT) 22, 23]. In particular, the frozen density 
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ular tool to treat intermolecular interactions in a simple 
and potentially exact embedding scheme. 

In the FDE approach the total electronic density of a 
system is partitioned into a set of subsystem densities 
(e.g., for two subsystems A and B: p = pA + Pb)- The 
total electronic energy is then written as 
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where Tg is the noninteracting kinetic energy functional, 
the external potential was partitioned as Uext = v^xt + 
vfy.^, and the superscript nadd denotes the nonaddi- 
tive energy contributions to the kinetic and exchange- 
correlation (XC) energies. The minimization of the en- 
ergy of Eq. (dJ with respect to one of the subsys- 
tem electron densities (e.g., pa), keeping the other one 
fixed and using the expansion of the electron density in 
terms of (auxiliary) Kohn-Sham orbitals tpf^, leads to the 
Kohn-Sham equations with constrained electron density 



The KSCED of Eq. (??) provide the solution for the 
ground state of subsystem A subject to the interaction 
with the frozen density of the subsystem B. If the 
KSCED of both subsystems are considered in a freeze- 
and-thaw procedure [25|, |37|, the full variational solu- 
tion for the total system, equivalent to the usual Kohn- 
Sham solution, is obtained, except for approximations 
eventually included in the nonadditive interaction terms. 
Henceforth, the acronym FDE wil be used to refer to this 
approach. 

So far, numerous applications of the FDE-based 
method, using the KSCED approach, have been pre- 
sented in literature to study the interaction of non- 
covalently bound (or weakly interacting) molecular com- 
plexes or solvated molecules |26l434| . In this cases in fact 
the nonadditive kinetic energy term can be efficiently ap- 
proximated by means of semilocal approximations |38l - 
41|. However, because in the conventional formulation 



the use of KSCED is limited to local/semilocal XC func- 
tionals, most studies have focused on hydrogen-bond and 
dipole-dipole systems, since these interactions are domi- 
nated by electrostatic and polarization effects that can be 
described relatively well at the semilocal level of theory 

IHEI. 

Recently the method has been also extended to the 
generalized Kohn-Sham (GKS) framework 4J, |45J, so 
that hybrid functionals can be considered as well. In 
this approach, named GKS-FDE, the nonadditive kinetic 
and XC terms in Eq. (|4]) are replaced by the nonadditive 
direct and residual interactions F^'^'^'^ and i?"*"^*^ which 
include respectively the orbital-dependent and the explic- 
itly density-dependent kinetic and XC terms. The pres- 
ence of nonlocal and/or orbital-dependent terms makes 
however the direct implementation of the method rather 
involved. Therefore, a practical computational scheme 
was proposed where these terms are approximated with 
local functionals 4J]. Thus, in practice in the GKS-FDE 
method the electronic structure is evaluated using a non- 
local exchange operator for the intra-subsystem XC ef- 
fects but with an embedding potential evaluated at the 
local/semilocal level of theory. A similar approximation 
was also proposed for the extention of the KSCED to 
the use of orbital-dependent optimized effective poten- 
tials 0. 

The GKS-FDE method was tested for a variety of hy- 
drogen bond and dipole-dipole interactions and proved 
to be a very good approach for embedded simulations 
with hybrid functionals 4J-|46|. In particular, with re- 
spect to FDE-based calculations employing semilocal XC 
functionals, it enables a general reduction of the embed- 
ding errors, because of the smaller overlap of subsystem 
densities at the hybrid level, and it allows a better de- 
scription of subsystem properties due to the removal of 
most of the self- interaction error (SIE) |46l448|. Further- 
more, it makes it possible to extend the application of 
FDE-based techniques to treat systems that are typically 
poorly described by semilocal XC functionals. 

An example of the latter is given by ground-state 
charge-transfer complexes (also called donor-acceptor 
complexes). These systems derive in fact a consider- 
able portion of their stabilization energy from an electron 
charge-transfer between two centers. Thus, unlike for hy- 
drogen bonds or dipole-dipole interactions, where electro- 
static and polarization effects dominate, the interaction 
in charge-transfer complexes is largely determined by or- 
bital interactions. As a consequence, the description of 
ground-state charge-transfer complexes is very challeng- 
ing for conventional semilocal XC functionals because of 
the SIE inherent in these methods. On the other hand, a 
good performance is generally obtained for hybrid func- 
tionals, especially when rather large portions of Hartree- 
Fock exchange are included [49|-|65[. Similar results also 
concern the charge transfer between organic molecules 
and metal substrates 66|, [67 [ . 



Due to the difScult application of DFT-based meth- 
ods using semilocal XC functionals to these systems, 
in the last years only few FDE investigations treated 
ground-state charge-transfer complexes and related prob- 
lems m, m [ei-iTi]. Some of them [4l|, [t^ made use 
of complicated or ad hoc reparametrized kinetic energy 
functionals in order to obtain improved kinetic embed- 
ding potentials and improve the description of the sys- 
tem. However, these kinetic functionals cannot be rou- 
tinely applied in FDE calculations for large systems. On 
the other hand, an extensive investigation of different 
kinetic energy functionals [68| showed that most popu- 
lar kinetic energy approximations perform similarly and 
reasonably well for complexes displaying small/moderate 
charge transfer, so that it appears reasonable to con- 
sider in such applications the same kinetic energy ap- 
proximation as those commonly used for hydrogen bond 
and dipole-dipole interactions. The role of the XC func- 
tionals in embedding calculations dealing with charge- 
transfer systems, instead, was not given much attention 
so far [36]. 

In this paper we cover this issue and demonstrate how 
hybrid functionals can be practically utilized within an 
embedding scheme to obtain an improvement of the per- 
formance in the description of embedded systems for a 
benchmark of charge-transfer complexes. 



METHOD AND COMPUTATIONAL DETAILS 



The KSCED 24, M and GKS-FDE HA 



ap- 



proaches are implemented in the TURBOMOLE pro- 



gram package [71|, version 6.4. All calculations were 
performed using the fZ^i? script, which uses a freeze-and- 
thaw procedure [37| to guarantee the full relaxation of the 
embedded ground-state electron density of both subsys- 
tems, until dipole moments of the embedded subsystems 
con verg ed to 10~^ au. A supermolecular def2-TZVPPD 



7% [73| basis set was employed in all calculations to ex- 
pand the subsystem electron densities. A monomolecular 
(or too small) basis set is in fact insufficient to carry on 
FDE calculations with good accuracy, resulting in basis- 
set errors larger than the FDE errors and the differences 
between the methods. The def2-TZVPPD basis set adds 
diffuse basis functions to the def2-TZVPP [73[ basis set, 
thus granting an accurate representation for the electron 
densities even at the relatively large bonding distances 
characteristic of the systems under consideration. Very 
accurate integration grids were employed to minimize nu- 
merical errors. Additional details about our implementa- 
tion and computational procedure are reported in Refs. 

ii-El 

In our investigation we considered a set of twelve repre- 
sentative ground-state charge-transfer complexes (C2H4- 
F2, NH3-F2, C2H2-CIF, HCN-CIF, NH3-CI2, H2O-CIF, 
NH3-CIF, HCN-NF3, HNC-NF3, HF-NF3, CIF-NF3, and 



C2H4-CI2). The geometries and the reference binding en- 
ergies were taken from Refs. |49l|7J, except for the C2H4- 
CI2 complex (for this latter case we optimized the geome- 
try at the MP2/aug-cc-pVTZ level of theory ^^^ and 
calculated the reference binding energy, in analogy with 



Ref. l7iasi?r^=i?F^^°(^) 



^^^MP2 



where E^, ^ ' 



is the binding energy computed at the CCSD(T)/aug-cc- 
pVQZ level of theory [zlH H^ and A^^P2 jg ^-^^ ^jf, 
ference between the MP2 binding energies computed at 
the complete basis set limit (56-extrapolation 81|) and 
with an aug-cc-pVQZ basis set; the resulting binding en- 
ergy is in very good agreement to that reported in Ref. 

The calculations were performed using conventional 
semilocal XC functional (BLYP [H, [^], PBE 84]) as 
well as hybrid methods (PBEO [H and BHLYP [82|,l83|). 
All calculations were performed using DFT-D3 dispersion 
corrections [86|- To compute the embedding nonadditive 
interactions, the nonadditive noninteracting kinetic en- 
ergy term was approximated using the GGA functional 
revAPBEk ^9, 40]. In the case of hybrid XC functionals 
also the nonadditive exchange term was approximated: 
the PBE and the BSSx [82| exchange functionals were 
employed for PBEO and BHYLP, respectively; the re- 
sulting methods are indicated as PBEO/PBE and BH- 
LYP/BLYP. In the case of semilocal XC functionals in- 
stead, no approximation for the nonadditive exchange 
contribution to the embedding potential is required: the 
same exchange functional is used for the embedding po- 
tential, subsystem and conventional supermolecule KS 
calculations. 

The quality of embedded densities was investigated by 
considering the absolute deviation of plane-averaged em- 
bedding density, defined as 



Ap{z) 



\pAix,y,z)+ pB{x,y,z)- p {x,y,z)\ dxdy 

(5) 

where p^iKS jg ^j^g electron density obtained from a (gen- 
eralized) KS calculation on the whole system. A quanti- 
tative measurement of the absolute error associated with 
a given embedded density was then obtained by comput- 
ing the embedding density error 



e 



1000 

N 



Ap{z)dz 



(6) 



with N being the number of electrons. In the evaluation 
of ^ only valence densities were considered to avoid nu- 
merical problems related to the high value assumed by 
the electron density in the core regions and in considera- 
tion of the fact that the core density is not important for 
the determination of chemical and physical properties of 
interaction between the subsystems, which are of interest 
here. 

The net charge on the i-th subsystem (monomer 
charge) was calculated for each method by considering 



a grid-based atoms-in-the-molecule partitioning 87| of 
the electron density into atomic basins flj as 



E 



/ p(r)rfr - Zj 



(7) 



with Zj the nuclear charge of the j-th atom and Ni the 
number of atoms in the i-th subsystem. For two inter- 
acting subsystems j^a + J^b = and thus we define the 
charge transfer as x = I^a| = Wb\- Reference values for 
the charge transfer were evaluated using relaxed MP2 
densities computed with an aug-cc-pVTZ basis set. 

Finally, to inspect the quality of the embedding ener- 
gies we considered the error defined as 

AE[pA,ps,p^^^] = i?'=-^pA,PB] - E^^^lp^^^] , (8) 

where E"'^^^ is the embedding energy, computed using the 
embedding densities into Eq. ([1]), and E'^^^ is the en- 
ergy resulting from a conventional (generalized) KS cal- 
culation on the full system. Further analysis was also 
performed considering, according to Ref. 45, the error 
decomposition 

AE[pA,PB,p'''^^] = Ar[p^,pB,P°^']+Ai?[p,pGKS]^ 

+AX5[p^,PB,P°^'], (9) 

where the error on embedding energy is partitioned into 
kinetic (AT), exchange {AX^), and relaxation (AD) 
contributions. For full details see Ref. Ua. 



RESULTS AND DISCUSSION 

In this section we report and analyze the results of 
the FDE-based calculations on the charge-transfer com- 
plexes. We evaluate the performance of hybrid and 
semilocal XC functionals (BHLYP, PBEO, BLYP, and 
PBE), elucidating what are the roles of the exact ex- 
change, of the geometry and of the charge transfer be- 
tween the subsystems. 



Embedding energy 

In Tab. |T]we report the binding energies obtained from 
different embedding methods (-E™*) and the correspond- 
ing results from conventional DFT calculations on the 
whole supermolecular system {Ej^ ). In addition, we 
report the deviation of these results from accurate refer- 
ence values (A™^*^ and A^^^, respectively) and the er- 
ror on embedding energies as defined in Eq. (??). We 
note that the embedding error on the total energy is ex- 
actly equivalent to the difference of the binding energies 
computed respectively from embedding and conventional 
supermolecular calculations (i.e. the DFT calculation 



on the whole system), since the contributions from iso- 
lated subsystems exactly cancel out. Moreover, all bind- 
ing energies reported in Tab. H] are dispersion-corrected 
and, because we use a supermolecular basis set, the basis 
set superposition error for the total-system calculation is 
treated in the same way in both DFT and FDE-based 
methods (Boys-Bernardi correction [88|). In Tab. Uwe 
also report in the last lines the mean signed error (MSE) , 
the mean absolute error (MAE), and the mean absolute 
relative error (MARE) with respect to the reference bind- 
ing energies. 

A first inspection of the results, focusing the atten- 
tion on E^^^ and A^^^, confirms the finding, well es- 
tabhshed in hterature [11, [HI, [13, IH [H, [H] , that hy- 



brid functionals outperform semilocal ones for the cal- 
culation of interaction energies of ground-state charge- 
transfer complexes. In fact, BHLYP with a MARE of 
16.63% is the most accurate method. However, when in- 
teraction energies from embedding calculations are con- 
sidered (i.e. E^"^^), all methods yield rather similar per- 
formances with MAREs in the range 30-37% (for a dis- 
cussion of the issues related to the comparison of embed- 
ding binding energies with experimental or theoretical 

HSl). This suggests that 



references see e.g. Refs. [26|, |2 



a strong error cancellation occurs for embedding calcula- 
tions using semilocal functionals, while the cancellation 
is much less in the case of hybrid approaches. The reason 
for this behavior can be traced back to the fact that, in 
general, conventional semilocal functionals overestimate 
significantly the binding energy of charge-transfer com- 
plexes 43. lia IsiI - Im . ISQl (MSE is largely positive) while, 
on the other hand, when used in methods based on the 
FDE theory they show usually a tendency to produce too 
low interaction energies (AE is negative in most cases). 
On the contrary, hybrid functionals are quite accurate 
for the calculation of the interaction energies of charge- 
transfer complexes and, at the same time, correspond- 
ing embedding calculations reproduce with good accu- 
racy the parent supermolecular GKS calculations. 

Concerning the error on embedding energies (Ai?), 
we note that the methods using hybrid XC functionals 
definitely outperform the ones employing semilocal XC 
approximations, yielding average deviations more than 
twice smaller. This improvement can be partially recon- 
duced, in analogy to what shown in Refs. |45|, |46|, to 
a smaller overlap of the subsystem densities (i.e. to a 
smaller SIE) at the hybrid level, that leads to a reduc- 
tion of the error on embedding energy. This effect is best 
recognized if we consider a set of systems where only one 
subsystem is varying, as for NH3 interacting with CI2, 
GIF, and F2 respectively. In this case the value of AE is 
increasing from 2.6 to about 7 mHa when the CI atoms 
are replaced by F atoms and this increase nicely corre- 
lates to the fact that the fluorine atom is affected by a 
larger self-interaction error than chlorine, so that it is 
plagued by a stronger overestimation of the diffuseness 



of the electron density at semilocal level of theory. How- 
ever, in general it is not possible to find a trend valid for 
all the systems due to the different details of the interac- 
tion when both subsystems are changed. 

In fact, despite the inclusion of exact exchange in the 
calculations may help to reduce the errors on embedding 
energies by lowering the effective overlap between the two 
subsystems, this is not the only effect determining the ac- 
curacy of the GKS-FDE energy calculations. In particu- 
lar for hybrid methods still a subtle error compensation 
between kinetic, relaxation, and exchange contributions 
occurs. Thus, to shed light on the different aspects de- 
termining the embedding energy error we can recur to 
the embedding energy decomposition introduced in Ref. 
45I . In Tab. |ll]we report the energy error contributions 
AT + AD and AX^ for all the hybrid methods (BHLYP 
and PBEO) and system considered. For each energy error 
contribution we report the MSE, MAE and MARE. The 
contributions due to kinetic (AT) and relaxation (AD) 
effects are reported summed together because both terms 
yield very large values (100-200% of the binding energies) 
and with opposite sign, so they only contribute to the to- 
tal error through a strong error cancellation. 

The analysis of the AT + AD contributions shows that 
the use of the subsystem electron densities from hybrid 
calculations can provide a significant improvement of the 
semilocal embedding energies (the MAEs are 1.35 and 
2.19 to be compared with 2.54 and 2.61), although the 
errors are still larger than the global ones reported in Tab. 
U (MAEs of 0.94 and 0.90 for BHLYP and PBEO, respec- 
tively) . This improvement is determined exclusively from 
the reduced overlap displayed by the densities obtained 
from hybrid calculations, since here no exchange approx- 
imation is considered. To make a quantitative evaluation 
of the effect it is possible to use the semilocal differential 
relative error (SDRE) statistical indicator, defined as 



N 

SDRE = — V 

N '^ 



1 A \AT, + ADA- lAEf^^^l 



E. 



■ref 



(10) 



with the sum extend over all the N systems reported in 
the table. A negative value of SDRE indicates that us- 
ing subsystem densities from a given hybrid calculation 
improves on average the calculation of kinetic and relax- 
ation energy contributions over the use of GGA densities. 
Using the SDRE indicator it is thus possible to compare 
energetic error terms with the same explicit level of ap- 
proximation, i.e. the semilocal kinetic energy approxi- 
mation and the density relaxation error [45| . Looking at 
the values of SDRE for the two hybrids considered in the 
table, it can be seen that in both cases these methods 
outperform the corresponding semilocal ones (negative 
values of SDRE). This result refiects the good quality of 
the electron densities computed with hybrid embedding 
methods (see next section) and evidences once more the 
importance of the density overlap in this context. In fact. 



TABLE I: Dispersion-corrected binding energies resulting from conventional supermolecular DFT (E^ 



and embedding 



{E^ ) calculations for several test charge-transfer complexes. The differences of the binding energies with respect to the 



reference ones are also reported (A 



Ei 



El"' and A^ 

praf 



E^, 



E^'' ) as well as the error on the embedding energy 



(AE; see Eq. (??)). The reference binding energies {E^"^ ) are reported in the second column. In all calculations we used the 
revAPBEk functional for the non-additive kinetic energies and a supermolecular def2-TZVPPD basis set. In the last lines we 
report the mean absolute error (MAE), the mean signed error (MSE), and the mean absolute relative error (MARE). All values 
are in mHa. 



System 



E7' 



e; 



■GKS 



Ei 



, GKS 



AE e: 



GKS 



. GKS 



AE 



BLYP 



BHLYP 



NF3-HCN 


1.67" 


0.96 


1.40 


-0.71 


-0.27 


0.44 


1.71 


1.11 


0.04 


-0.56 


-0.60 


C2H4-F2 


1.69'' 


5.21 


0.56 


3.52 


-1.13 


-4.65 


2.24 


0.11 


0.55 


-1.58 


-2.13 


NF3-HNC 


2.31" 


2.93 


3.07 


0.62 


0.76 


0.14 


2.56 


2.35 


0.25 


0.04 


-0.21 


C2H4-C12 


2.60" 


6.25 


4.76 


3.65 


2.16 


-1.49 


3.81 


3.91 


1.21 


1.31 


0.10 


NH3-F2 


2.88'' 


8.22 


1.16 


5.34 


-1.72 


-7.06 


3.51 


0.99 


0.63 


-1.89 


-2.52 


NF3-C1F 


2.92" 


4.14 


2.14 


1.22 


-0.78 


-2.00 


2.89 


1.40 


-0.03 


-1.52 


-1.49 


NF3-HF 


2.92" 


4.14 


3.22 


1.22 


0.30 


-0.92 


3.52 


2.71 


0.60 


-0.21 


-0.81 


C2H2-C1F 


6.07'' 


8.74 


5.09 


2.67 


-0.98 


-3.65 


6.96 


5.26 


0.89 


-0.81 


-1.70 


HCN-CIF 


7.74'' 


8.17 


6.70 


0.43 


-1.04 


-1.47 


7.90 


7.27 


0.16 


-0.47 


-0.63 


NH3-CI2 


7.78'' 


11.82 


9.23 


4.04 


1.45 


-2.59 


8.59 


8.67 


0.81 


0.89 


0.08 


H2O-CIF 


8.54'' 


10.30 


8.11 


1.76 


-0.43 


-2.19 


9.60 


8.77 


1.06 


0.23 


-0.83 


NH3-CIF 


16.92* 


25.13 


21.29 


8.21 


4.37 


-3.84 


21.03 


20.84 


4.11 


3.92 


-0.19 



MSE 


2.66 0.22 -2.44 


0.86 -0.05 -0.91 


MAE 


2.78 1.28 2.54 


0.86 1.12 0.94 


MARE 


71.47% 31.23% 72.56% 


16.63% 30.06% 32.43% 



PBE 



PBEO 



NF3-HCN 


1.67" 


1.42 


1.83 


-0.25 


0.16 


0.41 


1.48 


1.91 


-0.19 


0.24 


0.43 


C2H4-F2 


1.69'' 


5.39 


1.11 


3.70 


-0.58 


-4.28 


3.05 


1.13 


1.36 


-0.56 


-1.92 


NF3-HNC 


2.31" 


3.49 


3.61 


1.18 


1.30 


0.12 


2.97 


3.48 


0.66 


1.17 


0.51 


C2H4-C12 


2.60" 


7.61 


6.09 


5.06 


3.54 


-1.52 


5.96 


6.38 


3.41 


3.83 


0.42 


NH3-F2 


2.88'' 


8.73 


1.83 


5.85 


-1.05 


-6.90 


4.87 


1.89 


1.99 


-0.99 


-2.98 


NF3-C1F 


2.92" 


5.01 


2.87 


2.09 


-0.05 


-2.14 


3.59 


2.77 


0.67 


-0.15 


-0.82 


NF3-HF 


2.92" 


4.73 


3.82 


1.81 


0.90 


-0.91 


4.02 


3.80 


1.10 


0.88 


-0.22 


C2H2-C1F 


6.07'' 


10.56 


6.85 


4.49 


0.78 


-3.71 


8.84 


7.33 


2.77 


1.26 


-1.51 


HCN-CIF 


7.74'' 


9.63 


8.01 


1.89 


0.27 


-1.62 


8.46 


8.43 


0.72 


0.69 


-0.03 


NH3-CI2 


7.78'' 


13.48 


10.64 


5.70 


2.86 


-2.84 


10.90 


10.69 


3.12 


2.91 


-0.21 


H2O-CIF 


8.54'' 


11.63 


9.21 


3.09 


0.67 


-2.42 


10.21 


9.76 


1.67 


1.22 


-0.45 


NH3-CIF 


16.92' 


28.84 


24.40 


11.92 


7.48 


-4.44 


25.62 


24.31 


8.70 


7.39 


-1.31 



MSE 


3.88 1.36 -2.52 


2.17 1.49 -0.67 


MAE 


3.92 1.64 2.61 


2.20 1.77 0.90 


MARE 


91.53% 34.43% 71.65% 


45.82% 36.92% 29.83% 



" CCSD(T)-F12b/VTZ-F12, Ref. i4a. 

'' Wl benchmark calculations, Ref. [tJ. 

" CCSD(T) extrapolated to CBS limit, this work. 



the larger improvement is obtained for BHLYP (SDRE 
of -33.51%) due to its higher content of exact exchange 
with respect to PBEO (SDRE of -12.97%). 

Considering the exchange contributions to the embed- 
ding error {AX^) we note immediately that, for both 
methods, they are in general of opposite sign and of the 
same order of magnitude as the AT + AD terms. Hence, 



an important error cancellation effect occurs in all cases. 
This effect can be measured by considering the exchange 
differential error (XDE) and the exchange differential rel- 



TABLE II: Kinetic-relaxation (AT + AD), and exchange 
(AX ) contributions to the embedding energy error (see Ref. 
I45I for details) . In all calculations we used the revAPBEk func- 
tional for the non-additive kinetic energies and a supermolec- 
ular def2-TZVPPD basis set. At the bottom of each panel 
the mean signed error (MSE), mean absolute error (MAE) 
and the mean absolute relative error (MARE) are reported. 
For the AT-I- AD values also the semilocal differential relative 
error (SDRE; see text) is reported. For AX we computed 
the exchange differential error (XDE) and the exchange dif- 
ferential relative error (XDRE) [45l |. All values are in mHa. 



Systems 


BHLYP 


PBEO 




AT + AD 


NF3-HCN 


-0.4 


-0.4 


C2H4-F2 


2.6 


3.4 


NF3-HNC 


-0.3 


-0.2 


C2H4-CI2 


0.3 


1.1 


NH3-F2 


3.7 


5.5 


NF3-CIF 


1.1 


1.8 


NF3-HF 


0.6 


0.8 


C2H2-CIF 


2.2 


3.3 


HCN-CIF 


0.4 


1.2 


NH3-CI2 


0.9 


2.3 


H2O-CIF 


0.8 


1.9 


NH3-CIF 


2.9 


4.3 


MAE 


1.35 


2.19 


MSE 


1.23 


2.09 


MARE 


39.19% 


58.65% 


SDRE 


-33.51% 


-12.97% 




AX^ 


NF3-HCN 


1.0 


-0.0 


C2H4-F2 


-0.4 


-1.5 


NF3-HNC 


0.5 


-0.3 


C2H4-CI2 


-0.4 


-1.6 


NH3-F2 


-1.2 


-2.5 


NF3-CIF 


0.4 


-0.9 


NF3-HF 


0.2 


-0.6 


C2H2-CIF 


-0.5 


-1.7 


HCN-CIF 


0.3 


-1.2 


NH3-CI2 


-1.0 


-2.1 


H2O-CIF 


0.0 


-1.5 


NH3-CIF 


-2.7 


-3.0 


MAE 


0.72 


1.41 


MSE 


-0.32 


-1.41 


MARE 


19.10% 


34.10% 


XDE 


-0.41 


-1.28 


XDRE 


-6.79% 


-28.83% 



0.030 



0.025 



0.020 



'^ 



0.015 



0.010 



0.005 




«-oo9ro .8 .6 74 .2 0-2 
Distance (au) 

FIG. 1: Absolute deviation from the reference supermolec- 
ular electron density of plane-averaged embedding densities 
obtained from different embedding approaches, for the C2H4- 
F2 complex. The filled circles on the x-axis denote the atoms' 
positions. 



ative error (XDRE) statistical indicators [45 1 



XDE = 1^|A£;,|-|AT, + AA|, (H 
1=1 

1 A \AE,\ ^ \AT, + AD. 
M 2^ 



XDRE 



,ref 
'b,i 



(12) 



which provide (by increasingly negative values) an esti- 
mation of the improvement due to the error cancellation 
effect induced by the additional approximation on the ex- 
change nonadditive embedding terms. For BHLYP and 
PBEO the XDE and XDRE values are negative, indi- 
cating, as expected, an error cancellation between the 
kinetic-relaxation and the exchange contributions. How- 
ever, the XDE is three times larger for PBEO than for 
BHLYP, so in the former case a much more effective er- 
ror cancellation occurs. This explains why PBEO has the 
smaller energy error AE in Tab. HI despite it has kinetic- 
relaxation errors significantly larger than BHLYP. 



Embedding density 

In this section we consider the ability of different em- 
bedding methods to reproduce the electron density of 
supermolecular charge-transfer systems. This is an im- 
portant test for embedding approaches and provides di- 
rect insight into the quality of the embedding potential 
|. To discuss this issue we report 



[401,1411,144 

in Fig. [T]the plot of the absolute deviation of the plane- 
averaged densities (Ap; see Eq. (??)) for the C2H4-F2 
complex, as an example. For other systems similar plots 
are obtained (not reported). 



461,191 



The figure shows that, unUke for hydrogen bond sys- 
which are prototypical examples for FDE- 



44 



terns 

based studies, the density errors are quite different for dif- 
ferent methods, especially if hybrid and GGA approaches 
are compared. This trend can be well captured by con- 
sidering the integrated errors on embedding densities (^; 
see Eq. (??)), as resulting from different methods, which 
are reported in Tab. IIIII These data, along with the pro- 



TABLE III: Valence density errors ^ for GKS-FDE calcula- 
tions (BHLYP and PBEO) and FDE-based calculations using 
conventional BLYP and PBE functionals. In all calculations 
we used the revAPBEk functional for the nonadditive kinetic 
energies and a supermolecular def2-TZVPPD basis set. In 
the last line the mean absolute error (MAE) is reported. 



Systems 



Hybrid 
BHLYP PBEO 



Semilocal 
BLYP PBE 



NF3-HCN 


0.13 


0.24 


0.29 


0.29 


C2H4-F2 


1.58 


2.75 


6.94 


6.35 


NF3-HNC 


0.41 


0.49 


0.58 


0.58 


C2H4-C12 


3.34 


4.31 


5.83 


5.77 


NH3-F2 


2.24 


4.36 


9.84 


9.59 


NF3-C1F 


0.64 


0.99 


1.64 


1.73 


NF3-HF 


0.60 


0.74 


0.96 


0.94 


C2H2-C1F 


3.14 


4.32 


6.03 


6.02 


HCN-CIF 


1.76 


2.33 


3.13 


3.21 


NH3-CI2 


3.98 


5.48 


7.45 


7.60 


H2O-CIF 


2.37 


3.41 


4.89 


5.06 


NH3-CIF 


8.07 


9.37 


11.27 


11.15 



MAE 



2.35 



3.23 



4.90 



4.86 



portionality of the Ap(z) profiles reported in Fig. [T] (the 
embedding errors in the density have a similar spatial 
distribution) , show that for all systems the errors on em- 
bedded densities are directly correlated to the amount of 
exact exchange included in the calculations (see also next 
section). In fact, the best results are found for BHLYP, 
which includes 50% of Hartree-Fock exchange and has a 
MAE of only 2.35, comparable to that found in the case 
of hydrogen-bond and dipole-dipole interactions 44|, |46| . 
On the contrary, methods based on GGA XC functionals 
yield significantly larger errors, more than twice as large 
as the BHLYP ones. We note that such large variations of 
the embedding density error with the method are rather 
unusual and are not encountered in the case of hydrogen 
bond and dipole-dipole complexes 4Jj|46|. This indicates 
the special relevance of the inclusion of exact exchange 
for the treatment of charge- transfer complexes via FDE- 
based methods and fully agrees with the analysis made 
in the previous section about the role of relaxation effects 
for the energy errors of different methods. 

Beside the total quality of the embedding density, in 
the present context another quantity of great interest is 



TABLE IV: Charge transfer x resulting from conventional su- 
permolecular DFT and embedding calculations for several test 
charge-transfer complexes. The mean absolute error (MAE), 
the mean signed error (MSE), and the mean absolute relative 
error (MARE) respect to the reference MP2 charge transfer 
are also reported". 



Systems 



X 



X 



GKS cmb 

X A 



BLYP 



BHLYP 



MP2 



NF3-HCN 


0.017 


0.011 


0.012 


0.010 


0.009 


C2H4-F2 


0.082 


0.014 


0.029 


0.015 


0.018 


NF3-HNC 


0.007 


0.007 


0.006 


0.006 


0.005 


C2H4-C12 


0.093 


0.044 


0.061 


0.045 


0.053 


NH3-F2 


0.105 


0.018 


0.032 


0.015 


0.026 


NF3-C1F 


0.035 


0.013 


0.015 


0.009 


0.013 


NF3-HF 


0.025 


0.013 


0.018 


0.012 


0.018 


C2H2-C1F 


0.102 


0.051 


0.067 


0.049 


0.052 


HCN-CIF 


0.050 


0.028 


0.033 


0.025 


0.034 


NH3-CI2 


0.122 


0.060 


0.064 


0.044 


0.061 


H2O-CIF 


0.072 


0.034 


0.041 


0.028 


0.042 


NH3-CIF 


0.217 


0.209 


0.163 


0.126 


0.177 


MSE 


0.035 


-0.001 


0.003 


-0.011 




MAE 


0.035 


0.007 


0.005 


0.011 




MARE 


118.07% 


18.64% 


17.50% 


25.63% 






PBE 


PBEO 




NF3-HCN 


0.016 


0.011 


0.013 


0.011 


0.009 


C2H4-F2 


0.074 


0.013 


0.039 


0.013 


0.018 


NF3-HNC 


0.007 


0.007 


0.006 


0.006 


0.005 


C2H4-CI2 


0.093 


0.044 


0.074 


0.044 


0.053 


NH3-F2 


0.101 


0.017 


0.052 


0.015 


0.026 


NF3-CIF 


0.037 


0.013 


0.024 


0.011 


0.013 


NF3-HF 


0.025 


0.013 


0.020 


0.013 


0.018 


C2H2-CIF 


0.101 


0.050 


0.081 


0.049 


0.052 


HCN-CIF 


0.051 


0.028 


0.040 


0.027 


0.034 


NH3-CI2 


0.125 


0.061 


0.089 


0.050 


0.061 


H2O-CIF 


0.074 


0.035 


0.054 


0.030 


0.042 


NH3-CIF 


0.221 


0.206 


0.191 


0.162 


0.177 


MSE 


0.035 


-0.001 


0.015 


-0.006 




MAE 


0.035 


0.006 


0.015 


0.007 




MARE 


113.80% 


19.26% 


48.22% 


21.00% 





" The statistical indicators are computed using the error 
Ax*^ = X^ - x"^^^ with I =KS,GKS,emb. 



the net charge transfer predicted by different methods 
for the complexes under examination. The results for the 
various approaches, together with reference MP 2 values 
are reported in Tab. IIVI 

The table shows the well-known tendency of the 
GGA XC functionals to overestimate the charge transfer 
53,[5li, so that BLYP and PBE yield for x^^^ MAREs 
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well above 100%, while much better results (MARE of 
17.5%) are obtained by including a large fraction of ex- 
act exchange, as in BHLYP. However, a different situa- 
tion is observed when embedding calculations are consid- 
ered. These have, in fact, a marked tendency towards an 
underestimation of the charge transfer, compared to the 
corresponding supermolecular calculation on the whole 
system, probably related to the excessive repulsive char- 
acter of the approximated embedding potential. This 
behavior is also much more pronounced for GGA meth- 
ods than for hybrids. As a consequence, the BHLYP 
results are slightly worsened compared to the reference 
MP2 ones when embedding is used, while the semilocal 
methods display a strong improvement. 

In order to shine light on this finding, we consider that 
embedding densities computed at the semilocal level of 
theory, when used to compute the energy, provide higher 
energies for the complexes with respect to correspond- 
ing densities from supermolecular calculations (semilocal 
functionals stabilize the energy of charge-transfer com- 
plexes by overestimating the charge transfer to compen- 
sate for the absence of long-range interactions [49]). Be- 
cause, on the other hand, the isolated systems are not 
affected by this problem, this leads in general to a reduc- 
tion of the computed binding energy (see Tab. |l|, thus 
effectively compensating for the tendency of the semilocal 
XC functionals to over-bind. This explains the fact that 
in Tab. HI for semilocal functionals. A''™'' has a much 
smaller MAE and MARE than A'-'^^. It is worth noting 
that the same considerations also apply to the hybrid ap- 
proaches, but in these cases the effect is smaller. In fact, 
PBEO binding energies are slightly improved (on aver- 
age) when computed via an embedding approach, while 
BHLYP results end up to suffer from a slight underbind- 
ing. 



Role of the amount of exact exchange 

In previous sections we evidenced the important role 
of the exact exchange for the description of embedding 
energies and densities, by comparing the performance of 
FDE-based methods using semilocal XC functionals with 
that of GKS-FDE approaches using hybrid functionals 
based on exactly the same semilocal approximation. To 
push this investigation one step further here we consider 
more generally the family of one-parameter hybrid func- 
tionals 



E}^y^{a) = aEHF + (1 - a)E^'^^'' + E^ 



LYP 



(13) 



where Ehf is the Hartree-Fock exchange, a is a param- 
eter, and EVff is the Lee- Yang-Parr correlation energy 
functional [8J|. Eq. (??) reduces to BLYP for a = 0, 
to BHLYP for a = 0.5, and was already used to study 
the role of exact exchange in embedding calculations of 



hydrogen-bond and dipole-dipole interacting systems in 
Ref. S 

Fig. [2] reports, as a function of a, the errors on the 
embedding energies for selected systems as well as the 
individual contributions AT -I- AZ? and AX'^. For most 
of the systems the error /S.E is found to decrease with 
a until a ^ 0.5, and then to slightly increase for larger 
values of the parameter. This behavior is a consequence 
of an error cancellation between the AT -I- AD contri- 
bution (middle panel), which starts from a rather large 
positive value and decreases with a, and the AX^ term 
(right panel), which shows instead an evident parabolic 
shape with minima located around a — 0.3 and provides 
mostly negative contributions for a < 0.5. This behavior 
is quite different from the one observed for dipole-dipole 
and hydrogen-bond interactions (see Figure I of Ref. 45% 
which is generally characterized by positive contributions 
of AX , increasing with a. The present finding thus 
shows that for charge-transfer systems a more delicate 
balancing must be found between the need to increase 
the exact-exchange contribution, in order to reduce the 
density overlap, and the necessity to obtain the correct 
description of long-range XC interactions in the hybrid 
functional. 

The family of one-parameter hybrid functionals can be 
also be employed to study the effect of the exact exchange 
on the description of embedded densities. The result of 
this study was already partially anticipated in previous 
sections and it is clearly shown in Fig. [3l by increasing 
the amount of exact exchange the error on the embedding 
density is monotonically decreased. As mentioned in pre- 
vious discussions and in Ref. |46| this behavior is mainly 



due to the reduction of the overlap of the subsystem den- 
sities consequent to the reduction of the SIE. Of course, 
the effect is somehow attenuated for very large values of 
the parameter a because of the increased importance of 
the semilocal approximation used for the exchange term 
in the GKS-FDE scheme in these cases, as shown in right 
panel of Fig. [2l 

Finally, we can consider the ability of methods includ- 
ing different amounts of exact exchange to correctly re- 
produce the charge transfer. We thus report in Fig. |4]the 
difference Ax between the charge transfer computed for 
each value of a, with or without the embedding treat- 
ment, and the reference values obtained from MP2 re- 
laxed densities. The plot shows that considering full 
GKS calculations the charge transfer is generally over- 
estimated at the GGA level and better agreement with 
the reference is only found when a significative amount 
of Hartree-Fock exchange is considered. On the other 
hand, for embedding calculations the GKS charge trans- 
fer can be only reproduced when exact exchange is in- 
cluded, while calculations using a semilocal XC potential 
yield a strong underestimation with respect to the corre- 
sponding supermolecular calculation. In fact, the charge 
transfer computed from embedding methods appears to 
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FIG. 2: Errors on the embedding energy for selected systems are reported in the left panel as computed with the family of 
functionals of Eq. (??). The corresponding contributions AT-I- AD and AX are plotted in the middle and in the right panels, 
respectively. In all calculations we used the revAPBEk functional for the non-additive kinetic energies and a supermolecular 
def2-TZVPPD basis set. 
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FIG. 3: Errors on the embedding density for selected systems 
as computed with the family of functionals of Eq. (??). In all 
calculations we used the revAPBEk functional for the non- 
additive kinetic energies and a supermolecular def2-TZVPPD 
basis set. 



be almost independent from the percentage of exact ex- 
change included in the XC functional. To understand 
these results wc must consider that within FDE-based 
methods the number of electrons on each subsystem is 
fixed a priori, so that no real charge transfer between 
them can occur. The charge transfer must be instead 
mimicked by a strong polarization of the two subsys- 
tem densities, in such a way that their sum will give the 
same spatial distribution as the electron density of the 
full system. The role of the embedding potential is thus 
to provide a polarization of the subsystems, while each 
subsystem must respond properly to this perturbation. 
In the case of semilocal functionals they are well known 
to provide a poor description of polarization ^], thus 
a poor description of charge transfer must be expected. 
On the contrary, hybrid functionals have been shown to 
give a good response to polarization perturbations within 
the GKS-FDE scheme [4^ , thus they can provide a good 
reproduction of the charge transfer between the subsys- 
tems. 
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FIG. 4: Difference between the charge transfer computed with 
functionals defined by Eq. (??) and the reference MP2 value. 
In the left panel we report the results of GKS calculations 
on whole systems; in the right panel the corresponding em- 
bedding results are shown. In all calculations we used the re- 
vAPBEk functional for the non-additive kinetic energies and 
a supermolecular def2-TZVPPD basis set. 



Role of charge transfer 
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FIG. 5: Errors on embedding energy (AE), embedding den- 
sity (^), and the computed signed charge transfer on the di- 



halogen monomer GIF (i- 



,GKS 



) as a function of the 



charge transfer induced (at the GKS level) by the applied 
external field, for BHLYP and BLYP. The inset reports the 
external electric field used to generate different charge trans- 
fer values. 



In an analysis of the performance of embedding meth- 
ods for the simulation of ground-state charge-transfer 
complexes it is obviously interesting to study in details 
what is the relationship between the outcome of the cal- 
culations and the magnitude of the charge transfer. How- 
ever, this study cannot be performed straightforwardly 
by a simple comparison of different systems displaying 
different values of the charge transfer because in this way 
also many other details of the interaction are changed at 
the same time. For this reason in this section we prefer to 
focus the attention on a single system, namely the linear 
HCN-CIF complex, and force a modulation of the charge 
transfer through the application of an external constant 
electric field directed along the axis of the complex. 

In Fig. [5] we report the error on embedding energy 
(AE), the error on embedding density (^), and the error 
on the computed net charge on the dihalogen monomer 



GIF {i^' 



,einb 



,GKS 



) as a function of the charge transfer 



induced (at the GKS level) by the applied external field 
for the two functionals BHLYP and BLYP. The corre- 
sponding electric field is reported in the inset of Fig. [S] 
The figure shows that both functionals display the same 
behavior with the induced charge transfer, with all errors 
increasing as the charge transfer grows. In fact, the min- 
ima of the embedding errors are approximately located 
at zero charge transfer. The energy error decomposition 
(not reported) confirms that the increase of the energy 
error is essentially due to a large increase in the kinetic 
and relaxation contributions. 

Moreover, the lowest panel of Fig. [5] shows, in agree- 
ment with our analysis of previous sections, that the in- 



crease of charge transfer is in linear correspondence with 
an increase of the inaccuracy in its evaluation by embed- 
ding methods. This implies, as also shown by a compar- 
ison of the two upper panels of Fig. [5] with the lower 
one, that the embedding errors are correlated (almost 
linearly) with the inaccuracy of embedding approaches 
to correctly describe the amount of charge transfer. 



Role of geometry 

All previous investigations have been carried out at 
reference geometries. However, for charge-transfer com- 
plexes large differences in the equilibrium geometry can 
be expected for different methods. In particular, hybrid 
XG functionals are able to produce equilibrium structures 
very close to the reference ones but semilocal XG approx- 
imations lead to significantly shorter intermolecular dis- 
tances jSOl - lSa |54| . As shown in Ref . |4^ the geometry 
and the intermolecular distance play an important role in 
determining the accuracy of embedding calculations. It 
is therefore interesting, in the present context, to assess 
this issue for the limits represented by the structures op- 
timized at the semilocal BLYP and hybrid BHLYP level. 

To this end we consider in Fig. [6] the error on em- 
bedding energy (AE), the error on embedding density 
(^), and the error on the computed monomer charge 



(^' 



cinb 



,GKS 



) at different geometries interpolating lin- 



early between the BLYP and the BHLYP optimized ge- 
ometries of the complex G2H4-F2 (similar results are 
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FIG. 6: Errors on embedding energy (AE), embedding 
density (5), and the computed net charge on the ethylene 



monomer (v" 



,GKS 



) at different geometries interpolating 



between the BLYP and BHLYP equilibrium geometry of the 
C2H4-F2 complex. In the inset the embedding error decom- 
position for the BHLYP case is reported. 



found for other systems). The intermolecular distance 
between the F2 monomer and the plane of the ethylene 
is 1.91 A and 2.92 A for BLYP and BHLYP, respec- 
tively. The figure shows the important role of the ge- 
ometry for all the embedding results. This finding can 
be understood considering that at small intermolecular 
distances, i.e. at the BLYP geometry, the overlap be- 
tween the electron densities of the two subsystems rapidly 
grows. Moreover, because at reduced distances between 
the subsystems the orbital interactions are enhanced, also 
the charge transfer is highly increased, together with the 
error in its evaluation. As a consequence, large errors can 
be expected for the embedded densities, due to the well- 
known limitations of nonadditive kinetic potentials for 
strongly overlapping densities [23, |45|. Similarly, large 
errors are found for the embedding energies. In this case 
an energy decomposition of the BHLYP embedding error 
(see inset of Fig. ^ further confirms that the increase 
of the error is essentially due to a large increase of the 
relaxation-kinetic term, while only a minor role is played 
by the exchange term AX^ . Nevertheless, the latter be- 
ing always opposite in sign to the relaxation-kinetic con- 
tribution provides a moderate error cancellation which 
improves the final performance of the hybrid method. 

The results of Fig. [5] indicate in conclusion that, for 
the simulation of charge-transfer complexes via embed- 
ding approaches, the semilocal methods may display se- 
rious problems related to the use of relaxed structures 
(e.g. in geometry optimizations or molecular dynamics), 
beside the limitations already highlighted in the previ- 
ous sections. Therefore, the hybrid approaches appear 
definitely more robust in such applications. 



BHLYP/BLYP 




BLYP/BLYP 



Energy (Ha) 

FIG. 7: Density of states (DOS) for both GKS and hybrid BH- 
LYP (upper panel) and the semilocal BLYP (bottom panel) 
embedding calculations on the charge-transfer complex NH3- 
F2. The solid black line indicates GKS DOS of the whole sys- 
tem; the red dashed and solid turquoise lines indicate DOS 
from embedding calculations on the monomers NH3 and F2 
respectively. In the insets the GKS DOS (thick black line) is 
compared with the sum of the two subsystems DOSs (dashed 
magenta line). Energy values are in Ha. 



Orbital energies 

In previous sections we showed that only embedding 
approaches making use of hybrid XC functionals can re- 
produce the ground-state charge transfer of the investi- 
gated complexes with reasonable accuracy, while a strong 
underestimation occurs when semilocal approaches are 
used. This fact has important consequences on the accu- 
racy of embedding calculations, that nevertheless are par- 
tially, or even mostly, mitigated by error compensation 
effects. The improper description of the charge transfer 
can have however more subtle effects, especially on the 
one-electron spectrum of the subsystems. This is an im- 
portant issue, since the molecular orbitals of the embed- 
ded subsystems can be used as the basis for an analysis 
of the properties of the interacting subsystems and/or as 
input quantities in response calculations [ 96. . .97.] . 

The driving force for ground-state charge transfer in 
molecular complexes is a sizable difference in the chemi- 
cal potential between the constituting (isolated) subsys- 
tems, so that, when the subsystems start to interact, a 
fraction of the electron density must be moved from one 
fragment to another to equilibrate it. The final chemical 
potential of the complex is an average between those of 
the initial non-interacting subsystems and it is determin- 
ing the offset for the Kohn-Sham potential, thus de facto 
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fixing the energy of the molecular orbitals. In an ideal 
exact FDE-based embedding calculation the embedding 
potential shall be able, by definition, to yield for the em- 
bedded subsystems the same chemical potential as for the 
full complex and, of course, shall induce a polarization 
able to mimic the exact charge-transfer. As a side con- 
sequence, the single-particle spectrum of the subsystems 
will be also correct (as compared to the corresponding 
supermolecular calculation). Similar conclusions apply 
for any approximated embedding calculation having suf- 
ficient accuracy. On the other hand, if the embedding 
procedure fails to reproduce the correct charge transfer 
(i.e. correct with respect to the calculation on the full 
system), the polarization exerted by the embedding po- 
tential on the two subsystems will be insufficient, thus 
the chemical potential of the two fragments will be rather 
different from the one of the full system (too large in one 
subsystem and too low in the other). Consequently, the 
molecular orbitals of the two subsystems will lay at too 
high (for one subsystem) and too low energies (for the 
other). 

The two cases described above fit very well to the 
embedding calculations based on hybrid and semilocal 
functionals, respectively. This is nicely shown in Fig. [7] 
where, as an example, we report the density of states 
(DOS) of the NH3-F2 complex compared to that of the 
two embedded subsystems as resulting from BHLYP and 
BLYP embedding calculations. The comparison is in 
this case rather fair, since the two fragments compos- 
ing the NH3-F2 complex interact only weakly so that the 
molecular orbitals of the full system can be easily at- 
tributed to either of the constituent subsystems. The 
figure shows very clearly that BHLYP embedding calcu- 
lations yield subsystems orbitals which match well the 
ones from the full system (within 0.01 eV). Instead, the 
orbitals from the BLYP embedding calculations are sys- 
tematically shifted with respect to the ones of the cor- 
responding supermolecular calculation, with F2 orbitals 
shifted at lower energies and NII3 ones shifted at higher 
energies, consistently with the underestimation of charge 
transfer from NII3 to F2. Finally, we recall in favor of 
hybrid calculations that even supermolecular calculations 
using semilocal XC functionals provide a poor description 
of the single-particle spectrum, while better results can 



be only obtained by reducing the SIE [98|, |9 



CONCLUSIONS 

We applied methods based on the frozen density em- 
bedding theory to study the electronic properties of 
ground-state charge-transfer complexes and to clarify the 
ability of such approaches to treat accurately this class 
of systems and interactions. This topic is of particular 
interest due to the growing popularity of embedding ap- 
proaches in the field of computational chemistry and the 



importance of charge-transfer complexes in many chem- 
ical applications. However, to date, embedding calcu- 
lations have been mainly presented for hydrogen-bond 
and dipole-dipole interactions, while only few applica- 
tions dealt with systems where charge transfer plays a 
significant role [36|, |4l|, |68[ |69j . 

Our work clarifies some of the motivations for this lack 
of embedding calculations concerning charge-transfer sys- 
tems, showing that standard embedding approaches, 
making use of semilocal XC functionals, display several 
drawbacks in this context. Moreover, despite they can 
generally yield interaction energies of good quality (with 
respect to accurate references), this is due to a strong 
error cancellation effect and does not reflect a property 
of the corresponding Kohn-Sham supermolecular calcula- 
tion. For these reasons and because semilocal XC func- 
tionals are known to perform rather poorly for charge- 
transfer complexes even in conventional Kohn-Sham su- 
permolecular calculations, the use of such approaches 
does not appear appropriate. 

The introduction of the GKS-FDE method mitigates 
different problems in the FDE-based calculations and 
thus provides a valuable tool to handle embedding calcu- 
lations of charge-transfer systems. In fact, in the present 
work we showed the relevance of exact exchange contri- 
butions to reduce the embedding errors and improve the 
overall performance of the calculations, through a reduc- 
tion of the Coulomb self-interaction error and a better 
description of subsystems' electronic properties. 

In conclusion, FDE-based approaches making use of 
hybrid XC functionals appear as reliable and efficient 
methods to handle electronic structure calculations of 
systems characterized by charge-transfer interactions. 
This suggest the possibility of future studies in this di- 
rection making use of these tools to investigate charge- 
transfer interactions in complex environments. 
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